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The microscopic basis for the stability of itinerant ferro- 
magnetism in correlated electron systems is examined. To 
this end several routes to ferromagnetism are explored, using 
both rigorous methods valid in arbitrary spatial dimensions, 
as well as Quantum Monte Carlo investigations in the limit of 
infinite dimensions (dynamical mean-field theory). In partic- 
ular we discuss the qualitative and quantitative importance 
of (i) the direct Heisenberg exchange coupling, (ii) band de- 
generacy plus Hund's rule coupling, and (iii) a high spectral 
density near the band edges caused by an appropriate lattice 
structure and/or kinetic energy of the electrons. We furnish 
evidence of the stability of itinerant ferromagnetism in the 
pure Hubbard model for appropriate lattices at electronic den- 
sities not too close to half-filling and large enough U. Already 
a weak direct exchange interaction, as well as band degener- 
acy, is found to reduce the critical value of U above which 
ferromagnetism becomes stable considerably. Using similar 
numerical techniques the Hubbard model with an easy axis 
is studied to explain metamagnetism in strongly anisotropic 
antiferromagnets from a unifying microscopic point of view. 
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I. INTRODUCTION 

Even after several decades of research the microscopic 
foundations of itinerant ferromagnetism are not suffi- 
ciently understood. Indeed, in contrast to other col- 
lective electronic phenomena such as antiferromagnetism 
or conventional superconductivity there exists a remark- 
able gap between theory and experiment in this field. 
This has mainly to do with the fact that itinerant fer- 
romagnetism is a quantum- mechanical strong- coupling 
phenomenon whose explanation requires the application 
of non-perturbative techniques. Thus it belongs into 
the class of the most difficult many-body problems in 
condensed-matter physics. Significant progress was made 
in this field in the last few years due to the development 
and application of several new analytic and numerical ap- 
proaches. It is the purpose of this paper to present and 
discuss some of these new, exciting results. 

In Sec. 2 we recapitulate the derivation of a general 
lattice model for correlated electrons starting from a con- 
tinuum model, and discuss the various truncation steps 
which eventually lead to the Hubbard model. In partic- 
ular, the implications of the Lieb-Mattis theorem on the 
impossibility of ferromagnetism in these truncated mod- 
els in d = 1 dimension are analyzed. In Sec. 3 several 
microscopic mechanisms favoring ferromagnetism are dis- 
cussed. New and very recent results concerning ferromag- 



netism in Hubbard-type models obtained by various non- 
perturbative methods are presented and put into perspec- 
tive. Metamagnetic phase transitions and the need for 
the application of non-perturbative techniques for these 
investigations are the subject of Sec. 4. Finally, a con- 
clusion is presented in Sec. 5. 



II. ELECTRONIC CORRELATIONS AND 
MAGNETISM 

A. General lattice model 

Within the occupation number formalism the Hamil- 
tonian for electrons with spin a interacting via a spin- 
independent interaction V ee (r — r') in the presence_of 
an ionic lattice potential V mn (r) has the form 
H = H + H int , where 
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Here Vv(r), V>jr ( r ) are the usual field operators and 
n a (r) = ?/>+(r)?/v(r) is the local density. We note 
that the interaction term is diagonal in the space vari- 
ables r, r', i.e. it depends only on the (operator-valued) 
densities of the electrons at site r, r ' which interact 
via V ee (r — r'). The lattice potential entering the non- 
interacting part (|IJ) leads to a splitting of the parabolic 
dispersion into infinitely many bands which we enumer- 
ate by the index a. The non-interacting problem is then 
characterized by the Bloch wave functions Q k( r ) and the 
band energies e Q k- We may introduce Wannier functions 
localized at site Ri by 



ik Ri 



0akO) , 



(3) 



where L is the number of lattice sites, and thus con- 
struct creation and annihilation operators c^ ia ,c aia . for 
electrons with spin a in the band a at site Ri as 



ctia = J d3r Xai(r) V^(r) 
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Thereby the Hamiltonian may be written in lattice rep- 
resentation as 101 
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where the matrix elements are given by 
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We note that in contrast to the field-operator representa- 
tion defined in the continuum, the Wannier representa- 
tion does not lead to a site-diagonal form of the electron- 
electron interaction, i.e. the interaction does not only 
depend on the densities hi a = cf a c ia but contains ex- 
plicit off-diagonal contributions which will be discussed 
later. 



B. One-band models 

The Hamiltonian (||) is too general to be tractable in 
dimensions d > 1. Hence it has to be simplified using 
physically motivated truncations In particular, if the 
Fermi surface lies within a single conduction band, and if 
this band is well separated from the other bands and the 
interaction is not too strong, it may be justified to restrict 
the discussion to a single band (a = /3 = 7 = <5 = l). In 
this case (||) reduces to 



H l-band ^ ^ ^ij^itT^ 
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For most purposes this single-band Hamiltonian is still 
too complicated. Taking into account the weak overlap 
between neighboring orbitals in a tight-binding descrip- 
tion one may expect that the overlap between nearest- 
neighbors is most important. Hence the site-indices in 
(H) are restricted to nearest-neighbor positions. In the 
interaction this leaves us with a purely local contribu- 
tion Van — U, the Hubbard term, and the four nearest- 



neighbor contributions Vi 



V, Vi; 
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and Vujj = F' , which are o/f-diagonal in the site indices. 
The remaining one-band, nearest-neighbor Hamiltonian 
has the form || |f| 



H Uband - H Hub + V- 



NN 
1-band 



(9) 



H 



Hub 



-t 



E 

<2,i>>0" 



fiia&ja 



h.c.) + Uj2n l] h il (10) 



is the Hubbard model and 
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where 



is the contribution of nearest-neighbor interactions. Here 
hi = £<r^ and Si = V 2 E CTCT ' c^r CTCT ,c_,, where r de- 
notes the vector of Pauli matrices. In ([Tl]) the U-term 
describes a density-density interaction, the X-term is a 
bond-charge-site-charge interaction ( "density-dependent 
hopping"), the F'-term describes the hopping of local 
pairs consisting of an up and a down electron and, fi- 
nally, the .F-term corresponds to the direct Heisenberg 
exchange which is generically ferromagnetic in nature. 
The occupation number formalism illustrates particularly 
clearly that a spin-independent interaction plus the Pauli 
principle is able to lead to a mutual orientation of spins. 
Of all interactions in (jTJ) the Hubbard interaction U is 
certainly the strongest. Hence, in a final truncation step 
one may try to neglect even the nearest-neighbor inter- 
actions and retain only the on-site interaction U. This 
leaves us with the Hubbard model (fO), the simplest cor- 
relation model for lattice electrons 



1. The Hubbard model 

The Hubbard model was originally introduced in an 
attempt to understand itinerant ferromagnetism in 3d- 
transition metals H^jyJ . The expectation was that in 
this model ferromagnetism would arise naturally since 
in a polarized state the electrons do not interact at all. 
However, it soon became clear that in a ferromagnetic 
state the kinetic energy is also reduced. This makes 
the stability of ferromagnetism in the Hubbard model 
a particularly delicate problem. Indeed, the kinetic en- 
ergy with nearest-neighbor hopping usually favors anti- 
ferromagnetism. At half- filling (n — 1) and on bipar- 
tite lattices antiferromagnetism is a generic effect since 
it appears both at weak coupling (Hartree-Fock or Slater 
mean-field theory) and strong coupling (Anderson's "su- 
perexchange" mechanism). Hence it arises naturally in 
any perturbational approach and, in particular, is tract- 
able by renormalization group methods fl2]| . By con- 
trast, ferromagnetism is a non-trivial strong-coupling 
phenomenon which cannot be investigated by any stand- 
ard perturbation theory. 

The above discussion shows that, to understand the 
microscopic origin of itinerant ferromagnetism, non-per- 
turbative techniques are required. Unfortunately, there 
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are not many approaches of this type available; rig- 
orous mathematical methods (for recent reviews see 
ref. p3-15|), large-scale numerical methods &f^-ff9l, 
and variational approaches |20| - |25l are the ones most fre- 
quently used. 

As to rigorous results about ferromagnetism in the 
Hubbard model the Lieb-Mattis theorem [@ of 1962 is 
one of the most famous. It proves that for spin- and 
velocity-independent forces between electrons ferromag- 
netism cannot occur in one spatial dimension. This theo- 
rem applies to the general Hamiltonian (|5|) with infinitely 
many bands. On the other hand the general single-band 
model (|J) is not covered by the theorem unless (a) the in- 
teraction matrix element Vij mn is site-diagonal such that 
the interaction depends only on densities hi, and (b) the 
hopping and interaction does not extend beyond near- 
est neighbors. Hence the theorem applies to the model 
(|) with U,V ^0 but X, F, F' = 0, and thereby also to 
the Hubbard model (V = 0). We note that any critique 
of the single-band model (P) in view of the fact that it 
can lead to ferromagnetism in d = 1 in contrast to the 
Lieb-Mattis theorem would apply even more so to the 
Hubbard model since the latter is a particularly special 
single-band model. In other words: the fact that the 
Lieb-Mattis theorem applies to the Hubbard model but 
not to the general single-band model (@) does not make 
the Hubbard model a "more physical" model than (|^); 
after all it is only a special case of 

Another well-known theorem, that by Nagaoka |^(| of 
1966, provides explicit, albeit highly idealized, conditions 
under which ferromagnetism is stable in the Hubbard 
model with nearest-neighbor hopping. It proves that for 
U = oo the microscopic degeneracy of the ground state 
at half-filling (number of electrons N — number of lattice 
sites L) is lifted by a single hole, i.e. when one electron 
is removed (N = L — 1). In this case a saturated fer- 
romagnetic ground state is stable for any value of the 
hopping t on simple cubic and bec lattices, and for t < 
on fee and hep lattices. For the Nagaoka mechanism 
to work the lattice needs to contain loops along which 
the holes can move. Once the hole moves, the maximal 
overlap between the initial and final state clearly occurs 
in a ferromagnetic configuration. The problem is that 
Nagoaka's proof does not even extend to two holes, that 
a single hole is thermodynamically irrelevant, and that 
the limit of U = oo is highly unrealistic. 



III. MICROSCOPIC MECHANISMS FAVORING 
FERROMAGNETISM 

The Hubbard interaction is the result of an extreme 
truncation of the interaction in the general Hamiltonian 
(||). All interactions beyond the purely local part (e.g. 
nearest-neighbor density-density interactions, direct ex- 
change, band degeneracy and the associated Hund's rule 
couplings) are totally neglected. The Hubbard interac- 
tion is therefore very unspecific — it does not depend on 



the lattice at all and hence not on the spatial dimension. 
The lattice structure only enters via the kinetic energy. 
Therefore the stability of ferromagnetism in the Hub- 
bard model can be expected to depend in a sensitive way 
on the precise form of the kinetic energy |T^ , |l4| , p6| -p9| . 
Strategies to find the essential "kick" for ferromagnetism 
in the Hubbard model and more general models should 
then proceed in different directions: one may (i) keep in- 
teractions beyond the Hubbard-?/ (in particular the di- 
rect exchange term F in ([)])), (ii) keep band degeneracy 
and Hund's rule couplings, or (iii) find an appropriate 
kinetic energy and lattice structure. We will now discuss 
several recent results obtained along these lines. 



A. The importance of the Heisenberg exchange 
interaction 

The Heisenberg exchange interaction, caused by direct 
quantum-mechanical exchange of electrons at nearest- 
neighbor positions (the F-Term in ( pd| ) with F > 0), fa- 
vors the alignment of the electronic spins and hence sup- 
ports ferromagnetism in a straightforward way ]7|,p|,p^[ . 
However, since this interaction is rather weak (Hubbard 
estimated F ~ ^ eV for 3d-metals, such that F <C U) 
it cannot be the sole origin of itinerant ferromagnetism 
in systems like Fe, Co, Ni. Nevertheless it may be qual- 
itatively important, since it may well give a correlated 
system with more or less strong ferromagnetic tendencies 
the ultimate push and trigger ferromagnetism in spite of 
its smallness. It is therefore unjustified to neglect the ex- 
change interaction for merely quantitative reasons. This 
becomes particularly clear in the limit of large U (with 
U > |*|, |V|,iX|, F, F') close to n = 1, when the one- 
band model (0) can be transformed into an effective t-J- 
model |Q 
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where P projects onto the subspace without doubly oc- 
cupied sites. The effective exchange coupling 
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has an antiferromagnetic part, due to Anderson's su- 
perexchange but modified by the X-term, and a ferro- 
magnetic part, due to the direct Heisenberg exchange. 
Hence for large enough Heisenberg exchange F and/or 
Hubbard repulsion U the exchange becomes effectively 
ferromagnetic. This effect is completely neglected in 
the Hubbard model where even in the Nagaoka-limit 
(U = oo) the dimensionless parameter FU/t 2 is kept 
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1. Generalization of Nagaoka's theorem 

If the Heisenberg exchange coupling is taken into ac- 
count, it is possible to generalize Nagaoka's theorem to 
U < oo 0. We start with the Hamiltonian H?$ and , (|). 
This Hamiltonian, like the Hubbard Hamiltonian, com- 
mutes with the total spin S = J2i The eigenvalues 
of S 2 are denoted by S(S + 1). We are concerned with 
saturated ferromagnetic states with one hole below half- 
filling, i.e. with largest possible eigenvalue S max = N/2 = 
(L — l)/2. There are 2S max + 1 = L such states with 
the same energy eigenvalue. It can be shown that the 
ground states of Hnn with one hole (i. e. N — L — 1) have 
maximum total spin S = S max — (L — l)/2 and are non- 
degenerate (apart from the above-mentioned (2S rnax +l)- 
fold spin degeneracy) in the following cases |3l| ]: 

Case 1: On any lattice, if F > 0, t < and (a) X ^ t 

and U > C/ c (1) , or (b) X = t and U > C/ c (2) . 

Case 2: On lattices with loops, ifAT = t<0,F = 0, 
and U > Ui 2) . 

In both cases t > is allowed if the lattice is bipartite. 

These results are summarized in Table Bl The con- 

(1) (2) ^ 

stants Uc and Uc are given by 



U^ = Z(2\t\ + 

(x-tf 

F 

C/i 2) = z(2\t\ + 



V-F-2\t\ 

{X-tf 



F' — 



F 
F 



\F' 



(14) 
(15) 



where Z is the number of nearest neighbors. Hence, if 
F > ferromagnetic ground states are stable on any lat- 
tice for U larger than a finite critical value. For F — > + 
we have Uc — * oo, thus yielding Nagaoka's condition 
for the pure Hubbard model. This shows once more that 
the Heisenberg interaction F, which is neglected in the 
Hubbard model, provides an obvious mechanism for sta- 
bilizing ferromagnetic ground states at finite U. Note 
that since X and t are expected to be of the same order 
of magnitude, the sensitive dependence on F, due to the 
term (X — t) 2 /F, may cancel from . The dependence 
of U c on t, V, F is depicted in Fig. |l|. The case X = t is 
special, since in this case the stability of ferromagnetism 
can be achieved either by F > 0, or by F > and t < 
if the lattice has loops. Note that the case F > is 
not connected to the case F — by a limiting procedure, 
since only in the latter case the lattice is required to have 
loops. 

The critical couplings U^ and Uc are sums of terms, 
each of which corresponds to a typical energy scale. This 
means that the on-site interaction U has to be larger than 
the energy describing the paramagnetic state (bandwidth 
~ Z\t\), as well as the threshold energies for the onset 
of a charge-density wave or pha se separation (~ Z\V\), 
77-pairing superconductivity [p2| (~ Z|F"|), and a spin- 
density wave (~ (X — t) 2 /F). However, these terms do 



not enter separately, but appear in combinations, i. e. the 
effects interfere, as should be expected. 

The above conditions are sufficient conditions. The oc- 
currence of ground states with maximum spin outside the 
above parameter region is not ruled out. As in Nagaoka's 
theorem for the pure Hubbard model, the ferromagnetic 
ground state is an itinerant state with non-zero kinetic 
energy, but the proof of its stability cannot yet be ex- 
tended to doping beyond a single hole. 



2. Magnetic phase diagram within the dynamical mean-field 
theory 

Even for the simplest electronic correlation model, 
the Hubbard model, exact solutions are not available in 
d = 2, 3 dimensions, and numerical methods — whether 
exact diagonalizations or Quantum Monte Carlo (QMC) 
techniques — arc limited by the smallness of the systems 
that can be studied. Hence one would like to construct, 
at least, a thermodynamically consistent mean-field the- 
ory which is valid also at strong coupling. Such a (non- 
perturbative) approximation is provided by the exact so- 
lution of a model in d = oo. It is now known that in 
the limit d — * 00 |33| one obtains a dynamical mean- 
field theory for Hubbard-type models where the spatial 
dependencies become local, but all quantum fluctuations 
of the d-dimensional model are included p4]-|3^| . In fact 
the problem is equivalent to an Anderson impurity model 
complemented by a self-consistency condition p6| , ^7| , 
leading to dynamical mean-field equations which can be 
solved numerically, e.g. within a finite-temperature QMC 
scheme [^9); for reviews see refs. 40 4^]. We employed 
this numerical approach (for details see ref . jl2| ) to inves- 
tigate the influence of the direct exchange interaction F 
on the stabilization of ferromagnetism in the Hubbard 
model, using the one-band Hamiltonian (|9|)-([i~L]) with 
X = F' = 0j43). To take the limit d -> 00 the cou- 
plings in (|9|)-([llD have to be scaled appropriately p3| , |34| | , 
i.e. 



t 



F* 



v* 
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In the following we set t* = 1. In the limit d — > 00 
the V^*-term |3J] and F*-term reduce to their Hartree- 
contributions. Hence their influence is that of a gener- 
alized, i.e. spin- and site-dependent, chemical potential. 
In the homogeneous phase the spin- and site-dependent 
terms vanish. This implies that the nearest-neighbor 
interactions F* and V* become important only in the 
symmetry- broken phase. Consequently the susceptibili- 
ties of the model ([j|)-([ll]) with X = F' = may be calcu- 
lated from the pure Hubbard model [Q; this simplifies 
the matter considerably. 

The phase boundaries between the paramagnetic, an- 
tiferromagnetic and ferromagnetic phases may then be 
calculated in principle as follows [Q: (i) QMC simu- 
lations are performed in the homogeneous phase of the 
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pure Hubbard model for given U, temperature T, fill- 
ing n, and number of Matsubara frequencies A; (ii) for 
arbitrary values of F* , V* the appropriate susceptibil- 
ities are calculated; (iii) the inverse susceptibilities are 
extrapolated to A — * oo, if they are negative the ho- 
mogeneous phase is found to be unstable; (iv) to obtain 
ground state properties the calculated quantities have to 
be extrapolated to T — > 0. The results of these calcu- 
lations for n = 1, T = and a semi-elliptic DOS are 
collected in Fig. g, where the exchange coupling F* is 
plotted versus the Hubbard interaction U . We neglect 
the density-density terms in (fill), since they become im- 
portant only for V* - F*/2 > U @@. The solid line 
marks the phase boundary to the ferromagnetic phase. 
As expected the critical value of F* , F*, decreases with 
U: it depends on U as (F*(0) - F*(U)) cx U for small 
U (Hartree-Fock limit) and as F*{U) ocl/U for large U 
(Heisenberg limit). At U = 12 the value of F* neces- 
sary to induce ferromagnetism is seen to be as small as 
F* <~ 0.1. This shows how important even a weak ex- 
change coupling F is for the stability of ferromagnetism. 
For U > 4 there may well be a direct transition between 
the antiferromagnetic and the ferromagnetic phases. 

The shape of the DOS and the band filling n are very 
important factors concerning the stability of ferromag- 
netism, as will be discussed below (Sec. 3.3). Indeed, a 
symmetric DOS and filling n — 1, as used in the above 
calculation, are especially disadvantageous for ferromag- 
netism, partly because antiferromagnetism will be the 
generic magnetic order in this case. 



B. Band degeneracy and Hund's rule coupling 

Another important route to ferromagnetism may be 
taken by considering more than one energy band, namely 
by starting from M > 1 Wannier (or tight-binding) or- 
bitals. It is known from atomic magnetism that there are 
ferromagnetic couplings between electrons on the same 
atom leading to Hund's rules. These on-site "Hund's rule 
couplings" express the fact that by putting electrons in a 
maximum spin state an atomic exchange energy may be 
gained by the following mechanism. A spin wave function 
with maximum spin is always symmetric. Therefore, for 
the total wave function to be antisymmetric, the coordi- 
nate wave function must be antisymmetric. This reduces 
the probability for electrons to come close to each other 
which in turn lowers the Coulomb energy between them. 
This consideration establishes that in general there will 
be ferromagnetic on-site interactions even in a bulk sys- 
tem. Whether the presence of these terms is sufficient for 
ferromagnetism to appear in the ground state is, however, 
strongly dependent on their relative strength compared 
with the kinetic energy, on the lattice structure, electron 
density, etc. 

Let us therefore take a closer look at the terms in the 
Hamiltonian (g) in the case of M relevant bands. In this 



case we have to retain the band index a — I, ... M . Now 
there exist important on-site interactions even beyond 
the Hubbard interaction U = V^""", namely the follow- 
ing couplings that are off-diagonal in the band indices, 
and are hence only present for M > 1 bands: density- 
density interaction Vq 



^Hii i direct exchange inter- 



action F — , and hopping of double occupancies 

^0 = V?iif^- F° r simplicity we assume the orbitals to 
be equivalent, i.e. U is the same for all orbitals a, and 
Vo, Fa, Fq each have a fixed value for all pairs of orbitals 
a, (3. Furthermore it should be noted that for equivalent 
orbitals these parameters are not independent, but are re- 
lated by U = V + 2F , F = @. In addition to these 
Hund's rule couplings there are still the inter-site terms, 
namely the hopping tX which takes place only between 
like orbitals (this follows from the general derivation 
above), and the next-neighbor interactions V\ = Vy^" a , 

*i = V£g aa , Fi = v;;;;*-. f{ = v;;;;-. For simplicity, 

these next-neighbor couplings are assumed to be diago- 
nal in the band indices, i.e. they act only between like 
orbitals on neighboring sites. Finally, since we are deal- 
ing with equivalent orbitals, we assume that the next- 
neighbor parameters t, X\, F\, F[ each have the same 
value for all bands a. 

The resulting multi-band Hamiltonian then reads 
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is a straightforward generalization of the one-band Ha- 
miltonian (^|) to more than one band, and 
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i;a<(3 

1 
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+ F o(c»TC 4 + aiC l(3i c l( g T + h.c.) 



(19) 



Several of the processes contained in the Hamiltonian 
( p7| ) are illustrated in Fig. ||. 

The physical picture of bulk ferromagnetism, originally 
put forward by Slater, is the following. If the on-site 
Hund's rule couplings are strong enough, they lead to 
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an independent ferromagnetic alignment of spins on each 
atom. In this situation the kinetic energy plays a de- 
cisive role since it can serve to communicate the spin 
alignment across the solid. Indeed, if the alignment on 
neighboring atoms were different, the hopping of elec- 
trons from one atom to the next would generate on-site 
Hund's rule interactions and thus increase the energy. 
These interactions can only be avoided if the spin align- 
ment on neighboring atoms is the same, implying global 
ferromagnetism. This mechanism for the stabilization 
of ferromagnetic order works all the better the larger 
the number M of orbitals, i. e. bands, is onto which 
an electron can hop; it does not work for a single band 
(M = 1). Therefore one should expect band degeneracy 
to favor ferromagnetism. This very qualitative picture 
is indeed found in Stoner mean- field theory J46| |: Since 
the non-interacting DOS is proportional to the number 
of degenerate bands, the critical interaction decreases as 
U^ toner (M) ~ 1/M. On the other hand, few exact re- 
sults are known in the case of degenerate bands. For 
example, for a one-dimensional chain with two orbitals 
(M = 2) and infinite on-site Coulomb interactions it has 
been shown that the ground state is ferromagnetic for 
N = L + 1 electrons p7| . This statement has been ex- 
tended to L + 1 < N < 2L — 1, i.e. up to one electron 
less than half-filling The important point is that 

on a one-dimensional chain with only one orbital and 
hopping between nearest neighbors, Nagaoka's theorem 
is not applicable, since the lattice does not have loops. 
In that case the ground state is degenerate with respect 
to the total spin S. On the other hand, if there are two 
orbitals the loop property is restored, and the ferromag- 
netic states with maximum spin become the only ground 
states. 

Furthermore, the following rigorous result can be es- 
tablished for the Hamiltonian (17) at half- filling |l9| . For 
TV = ML electrons, the ground states of Flu-band have 



maximum spin S = S„ 



ML/2 if 



2V >F > 



U, 



(1,2) 



l + M/2 



and Fi > 



(20) 



(12) 

where t/ c ' are the critical values for U in the single 
band system, as given in (|l4|)-(|l5|). The meaning of these 
bounds is the following. The requirement that 2Vq > Fq 
and F be greater than a certain threshold leads to an 
alignment of the electronic spins on an isolated atom. 
On the other hand, ferromagnetism within each band is 
brought about by the next-neighbor exchange F\ > and 
the Hubbard interaction U larger than a threshold related 

(1 2) 

to Uc ' ■ The combination of these two effects (using 
the fact that U = Vb + 2Fo for equivalent bands) leads 
to a critical value for the Hund's rule coupling F , which 
indeed becomes lower as the number of bands increases. 

While this result contains some ideas of Slater's pic- 
ture, it does not explain the itinerant aspects of multi- 
band ferromagnetism, since at half-filling the ferromag- 
netic ground states are insulating. So far, this result can 



only be modified to apply also to Nagaoka's case 
hole, N — ML — 1) §§. 



C. Kinetic energy and lattice structure 

The stability of ferromagnetism is intimately linked 
with the structure of the underlying lattice and the 
kinetic energy (i.e. the hopping) of the electrons 
fp7[^§6}j2|]. This is supported by several facts: (a) Na- 
gaoka's proof of ferromagnetism in the Hubbard model 
for a single hole at U — oo |2(| depends on the existence 
of closed loops along which the hole can move, and (b) on 
bipartite lattices emizferromagnetism is the generic mag- 
netic state making it hard for ferromagnetism to become 
stable. Hence non-bipartite lattices with loops (or with 
a kinetic energy involving hopping between nearest and 
next-nearest neighbors sites effectively leading to a mo- 
tion on loops) should be expected to support ferromag- 
netism because the competing antiferromagnetic tenden- 
cies are severely weakened, and because the correspond- 
ing DOS of non-interacting electrons is asymmetric and 
thus has a peak at a non-symmetric position. Indeed, a 
peak at one of the band edges as in the case of the fee 
lattice is favorable for ferromagnetism |l0|]l6|Jl7| , p0| p2[ . 
This is supported by the observation p2| that Co and 
Ni, having non-bipartite hep and fee lattice structure, 
respectively, show a full magnetization while bcc-Fe has 
only a partial magnetization. 



1. A model density of states 

To gain insight into why a DOS with a peak (or more 
precisely with a high spectral density) at the band edge, 
e.g. at the lower band edge for n < 1, may be favorable 
for ferromagnetism, we study non-interacting electrons 
with the following model DOS: 



(21) 



where A is the width of the band and A parameterizes 
the asymmetry. The lower band edge is — (A + 3) A/6 
such that the first moment of N°(E) vanishes. For A — 
— 1(+1) the DOS has a triangular shape with the peak 
at the lower (upper) edge, while for A — it is flat. We 
wish to calculate the energy difference Se between the 
fully polarized and the paramagnetic state 



5e = ef erro (n,A) 



i(n,A) 



(22) 



as a function of A and the band filling. It is easy to con- 
firm that for all n 5e is lowest for A = — 1. The reason is 
this: In the paramagnetic state N non-interacting elec- 
trons (N/2 electrons with spin up and down each) fill the 
lowest N/2 /c-states up to an energy E^ ra , while in the 
ferromagnetic state the N singly occupied fc-states below 



G 



E ferro > ^para ^ j g ^ ^ £_ states with en _ 

ergy above _E^, a ™ are also occupied. The higher the DOS 
is at the lower edge the less the additional N/2 fe-states 
are forced into high-energy states, i.e. the lower Ep Tro 
will be. Then 5e is kept at a minimum. 



2. The Hubbard model on fee type lattices 

As explained in the previous section, a lattice structure 
which gives a high DOS at low energies is expected to 
be favorable for ferromagnetism (at n < 1) because this 
situation reduces the loss in kinetic energy. The question 
remains: is a strongly peaked DOS sufficient to induce 
ferromagnetic order in the single band Hubbard model 
without additional interactions? 

In this section we discuss the results of a QMC-investi- 
gation of the stability of ferromagnetism in the single 
band Hubbard model in the limit of infinite dimensions. 

Based on the considerations in Sec. 3.3.1, the fee lattice 
is expected to be a good candidate for ferromagnetism 
because of its high (divergent) DOS at the lower band 
edge. The fcc-lattice can be g eneralized to higher dimen- 
sions in different ways [p^pTf . Here we use the definition 
of an fee lattice as the set of all points with integer cubic 
coordinates summing up to an even integer [p0|] . It is 
a non-bipartite Bravais lattice for any dimension d > 2. 
For d = 2 it is identical to the square lattice. Nearest 
neighbors are connected by two different unit vectors on 
a simple hypercubic (he) lattice. The coordination num- 
ber is hence Z — 2d(d — 1). With the proper scaling of 
the hopping term, (Jlq), the non-interacting DOS of the 
generalized fee lattice can be calculated in d — oo p3| as: 



N°„, 



(E) = e-( 1+ ^ E V 2 /yJn(l + y/2E) 



(23) 



which has a strong square-root divergency at the lower 
band edge, — 1/^/2, and no upper band edge. Since we 
choose a positive hopping integral {—t > 0), N^ CC (E) 
might be regarded as the DOS of holes rather than elec- 
trons. A full polarization of holes would hence corre- 
spond to a maximal (though not full) polarization in a 
more than half-filled band. 

While the three dimensional fee lattice has no square- 
root but only a logarithmic divergency at the band edge it 
is worth mentioning that a square-root divergency arises 
on the fee lattice in any dimension if there is an addi- 
tional next nearest neighbor hopping of the size t' = t/2 
between sites that are linked by two unit vectors in the 
same direction on the he lattice. The energy dispersion 
and the DOS of this so-called "half-hypercubic" (hh) lat- 
tice [^lj can easily be obtained from the he lattice as: 



£hh(k) 



KM 



2t 



h c 



-4M 



3t 



K c N° hc (VE-- 



■3t) 



VE + 'St 



(24) 
(25) 



In the limit d — > oo, the hh and fee lattices become equiv- 
alent. 

The fact that the fee lattice provides a good 'envi- 
ronment' for ferromagnetism has also been supported by 
variational studies of the stability of the Nagaoka state 
po| , p2| . Variational calculations provide limits for the 
critical density n c and the critical interaction U c where 
saturated ferromagnetism becomes unstable: While for 
the he lattice in d = oo the stability regime shrinks to 
the point U c = oo, S e = @, for the fee DOS (§|) there 
is a critical line U c (n) with U c (0) = and U c (l) = oo 
HI . The Nagaoka state is always unstable in the case of 
electron doping (n > 1). 

Recently, these variational boundaries were qualita- 
tively confirmed by Uhrig's calculation of the exact sin- 
gle spin-flip energy of the Nagaoka state in d — oo . 
While on the hh lattice U c vanishes at low densities, U c 
remains finite for all densities in the case of the "lami- 
nated" lattice which is a different generalization of the 
fee lattice without a divergent DOS at the band edge. 
His results, too, emphasize the subtle dependence of the 
stability of ferromagnetism on the lattice structure. 

Antiferromagnetism is not expected on the fcc-lattice 
in high dimensions because the difference of the numbers 
of not frustrated bonds and frustrated bonds is only of 
the order of d resulting in an effective field of the order 
of t 2 d oc 1/d J52|. Even in d = 3 antiferromagnetism is 
frustrated and is expected only very close to half-filling. 

To detect a ferromagnetic instability we calculated the 
temperature dependence of the uniform static suscep- 
tibility, \F, from the two-particle correlation functions 
|f42f . At an intermediate interaction strength of U = 4 
we found the ferromagnetic response to be strongest 
around quarter filling (n ~ 0.5). xf obeys a Curie- 
Weiss law (Fig. Q) and the Curie temperature T c can 
safely be extrapolated from the zero of x^ 1 to a value 
of T c = 0.051(2) at n = 0.58. Below T c the magneti- 
zation m grows rapidly, reaching more than 80% of the 
fully polarized value (m max = n = 0.58) at the lowest 
temperature which is only 30% below T c . The three data 
points m{T) (Fig. |) are consistent with a Brillouin func- 
tion with the same critical temperature of T c = 0.05 and 
an extrapolated full polarization at T — 0. A saturated 
ground state magnetization is also consistent with the 
single spin-flip energy of the fully polarized state which 
is positive at the present parameter values |^l| . 

Translated to the three dimensional fee lattice with 
Z = 12 nearest neighbors and a bandwidth of W = 16t, 
the critical temperature becomes T C (3D) « 0.011IV fl53| |. 
Thus, despite the oversimplifications of the single band 
Hubbard model, the resulting Curie temperature has a 
realistic order of magnitude of 500-800K for typical values 
of W around 5 eV. 

In order to answer the question if the system is still 
metallic we also calculated the single particle spectrum 
f34j . We find that the system is metallic since both spec- 
tra have a finite value at the Fermi level. 

We conclude that within the dynamical mean-field the- 
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ory a DOS with sufficiently large spectral weight at low 
energies is able to induce itinerant ferromagnetism in the 
single band Hubbard model. Detailed investigations of 
the dependence of ferromagnetism on the electronic den- 
sity, the interaction strength, and in particular the exten- 
sion to more realistic densities of states (e.g. fee lattice in 
d — 3) is under progress and will be presented elsewhere 
& 

In the d = oo limit the dynamics of the correlated 
system is fully taken into account, while spatial fluc- 
tuations are suppressed p4|-|3q]. One might therefore 
suspect that the stability of ferromagnetism is somehow 
overestimated in this approach (in particular since the 
rivaling antiferromagnetism is completely absent on the 
fee lattice). However, most recently similar results were 
reported for the Hubbard model with nearest and next- 
nearest neighbor hopping t and t', respectively, in dimen- 
sions d = 1 and d = 2 E3| , which are consistent with 
the results in d = oo. Hence the existence of itinerant fer- 
romagnetism in the pure Hubbard model with a suitable 
kinetic energy seems to be established at last. 

IV. METAMAGNETIC PHASE TRANSITIONS 

An issue related to the stability of (itinerant) ferromag- 
netism in the Hubbard model is the question concerning 
metamagnetism |35| in this and other models. Here it 
is an external magnetic field H which helps to suppress 
the (not necessary long-range) antiferromagnetic corre- 
lations in the system and thereby induces a pronounced 
transition from a state with low magnetization to one 
with high magnetization. At the metamagnetic transi- 
tion the magnetization curve m{H) shows an up-turn 
such that the susceptibility x{H) — dm/dH has some 
kind of maximum. This feature serves as a convenient 
general definition for "metamagnetism" . Metamagnetic 
transitions were first observed in strongly anisotropic 
antiferromagnets of which FeCl2 and DyaAlsO]^ (DAG) 
are well-studied prototypes ]56]]. These materials are in- 
sulators where the valence electrons are localized at the 
Fe and Dy ions, respectively. The arising local moments 
order antiferromagnetically and are strongly anisotropic 
in the sense that they are constrained to lie along an 
easy axis e. In this case a spin-flop transition in an ex- 
ternal magnetic field H || e cannot occur. Apart from the 
above materials there are also conducting systems that 
most probably belong to this class, e.g. the conductors 
XJAi- x B x (where A = P, As; B = S, Se) g], SmMn 2 Ge 2 
|3 and TbRha-xIr^Sia §|. 

Hitherto completely different theories are employed to 
describe metamagnetic phase transitions in these differ- 
ent, strongly anisotropic antiferromagnets. Investiga- 
tions of localized systems are usually based on the Ising 
model, where more than one interaction has to be intro- 
duced to describe the experimentally observed first or- 
der phase transitions f39|. With antiferromagnetic cou- 



pling J between the Z nearest-neighbor (NN) spins and 
a ferromagnetic coupling J' between the Z' next-nearest- 
neighbors (NNN) one obtains 

Hising = J SiSj — J' SiSj — 2H Si . (26) 

NN NNN i 

In Weiss mean-field theory this model shows two different 
types of phase diagrams depending on the parameter R = 
Z'J'/{ZJ) @. For R > 3/5 the first and second order 
phase transition line join smoothly at the same point, 
producing a tricritical point (TCP), while for R < 3/5 
there is no common endpoint (see Fig. ^|). However the 
scenario of Fig. ^3 was not found when evaluating ( |2"6| ) 
beyond mean-field theory . 

For itinerant electron metamagnetism (IEM) Moriya 
and Usami J62| proposed a Landau theory, where the 
parameters have to be deduced from microscopic models. 
The idea is to calculate first the independent electron 
band structure and to introduce the Coulomb interaction 
within the random phase approximation. 

It is our purpose to investigate the origin of metam- 
agnetism in strongly anisotropic antiferromagnets from 
a microscopic, quantum-mechanical point of view, and 
to describe different kinds of metamagnets (i.e. metallic 
and insulating, band-like and localized systems) qualita- 
tively within a single model. To this end we study the 
Hubbard model (110) with the additional constraint that 
the antiferromagnetic magnetization m st lies parallel to 
the external magnetic field H |6^]. In this way the exis- 
tence of an easy axis e along which H is directed, such 
that e || m st || H, is incorporated in a natural way. By 
this approach, both kinetic energy and Coulomb interac- 
tion are captured microscopically, whereas the relativistic 
corrections (responsible for the easy axis) are not. This 
procedure is justified since the relativistic corrections are 
of 0(1O~ 2 eV) and are thus small compared to the kinetic 
and Coulomb energy which are of O(leV). Therefore the 
existence of an anisotropy axis e and the correlations de- 
scribed by the Hubbard model are quite unrelated. Note, 
that the existing Ising and IEM theories do not treat the 
kinetic energy and Coulomb interaction microscopically, 
but within an effective model. 

A perturbative treatment of the Hubbard model with 
easy axis in the weak and strong coupling limit shows 
that the appearance of a tri- or multicritical point is 
a delicate matter, since neither of these two limits is 
able to describe a change of the transition from first to 
second order j64|. Apparently the entire transition sce- 
nario depends sensitively on the value of the electronic 
on-site interaction U. To study this point in greater 
detail we have to go to intermediate coupling. In this 
non-perturbative regime we employ again QMC simula- 
tions to calculate the magnetization m(H) and the stag- 
gered magnetization m st (H) of the Hubbard model ( J10| ) 
in d = oo pi] . As the results do not much depend 
on the precise form of the density of states we choose 
N°{e) = [(2t*) 2 - e 2 ] 1/2 /(27rf* 2 ), setting t* = 1 in the 
following. All calculations are performed at half-filling. 
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The results for m(H) and m s t (H) are used to construct 
the H—T phase diagram at U — 4 (Fig. ||). It displays all 
the features of Fig. ||b. In particular, the first order line 
continues into the ordered phase, separating two different 
AF phases: AFj (where m ~ 0) and AF77 (where m > 
0). The position of its endpoint cannot, at present, be 
determined accurately (dotted line) . This phase diagram 
is surprisingly similar to the experimental phase diagram 
of FeBr 2 f§. 

A change in U and the filling n will affect the phase 
diagram quantitatively and qualitatively. These results 
will be reported elsewhere [p6[. 



V. CONCLUSION 

In the last few years, and especially most recently, con- 
siderable progress was made in our understanding of the 
microscopic origin of itinerant ferromagnetism. These re- 
sults were obtained on the basis of well-defined lattice 
models of correlated electrons, of which the one-band 
Hubbard model is a particularly important ingredient, 
by applying new, non-perturbative techniques, ranging 
from rigorous to large-scale numerical methods. There 
exists convincing evidence now that on appropriate, non- 
artificial lattices, or for an appropriate kinetic energy, 
itinerant ferromagnetism is stable even in the pure Hub- 
bard model, for electronic densities not too close to half- 
filling and large enough U . Important ingredients are: 

(i) lattices with loops (or a kinetic energy allowing for 
motion on loops, e.g. with t, t' hopping) such that the Na- 
gaoka mechanism works and antiferromagnetism is sup- 
pressed, 

(ii) a large spectral weight near the band edge. (We note 
that this condition goes far beyond the mean-field Stoner 
criterion for ferromagnetism, UN(0) — 1, where N(Q) is 
the DOS at the Fermi energy.) 

The direct exchange interaction, as well as band degen- 
eracy, will strongly reduce the critical value of U above 
which ferromagnetism becomes stable. 

Furthermore, using the dynamical mean-field theory to 
solve the Hubbard model with easy-axis in the intermedi- 
ate coupling regime, it appears to be possible to describe 
various, different forms of metamagnetism within a single 
microscopic theory. 

The development and application of controlled, non- 
perturbative techniques will continue to be of particular 
importance for the investigation of correlated electron 
systems. 
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TABLE I. Sufficient conditions for ferromagnetic ground 
states with one hoie. 



. Y rr . rr (i) any t bipartite lattice 

la t > U, A 7^ t, U > Uc i - n u • , ... 

t < non-bip. lattice 

,, „ n v 4- t t rr(2) any t bipartite lattice 

lb p > U, A = t, U > U c „ ,. . . 

f < non-bip. lattice 

p n v + rr ^ rr< 2 ) * ^ bipartite, with loops 

1 t — u, a = t, ty > u c n , . . , . 

f < non-bip., with loops 



FIG. 1. Critical value U c vs. exchange interaction F for 
different V, X, and F' = 0. For U > U c the ground state 
is ferromagnetic. See Table | for details. 

FIG. 2. Critical value of the direct exchange coupling F* 
above which ferromagnetism is stable vs. Hubbard interac- 
tion U. The phase boundaries were calculated from the diver- 
gence of the ferromagnetic (diamonds) and antiferromagnetic 
(squares) susceptibilities, respectively. The QMC results were 
extrapolated to T = 0, the filling is n = 1 (F: ferromagnetic 
phase, AF: antiferromagnetic phase, P: paramagnetic phase). 
Dotted line: Hartree(-Fock) theory; dashed line: Heisenberg 
limit. 

FIG. 3. Illustration of the interaction and kinetic energy 
in a two-band model: U, Vq, Fq are the Hubbard interaction, 
the density-density interaction, and the direct exchange inter- 
action, respectively, acting on the same site, while Vi and Fi 
act between neighboring sites i and j, and t is the hopping. 

FIG. 4. Magnetization m (diamonds) and inverse ferro- 
magnetic susceptibility Xf- 1 (circles; values multiplied by a 
factor of 4 to use the same scale) for U — 4 and n = 0.58. 
Error-bars are of the size of the symbols or smaller. (Note 
that the value of Xf 1 at T = 0.05 is a data point, not an ex- 
trapolation.) The dotted line is a linear fit to Xp 1 , the dashed 
line a fit with a Brillouin function to m. 

FIG. 5. Schematic H — T phase diagram for a) a typical 
Ising-type metamagnet (TCP: tricritical point), b) the Ising 
model (2^) in mean-field theory with R < 3/5 (CE: criti- 
cal endpoint, BCE: bicritical endpoint) p3]. Full lines: first 
order transition, broken lines: second order transition; AF: 
antiferromagnetic phase, P: paramagnetic phase. 

FIG. 6. H — T phase diagram for the d — oo Hubbard 
model with easy axis along H at n = 1 and U = 4 |j4|; same 
notation as in Fig. Hb. 
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